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Abstract 



We show that current estimates of the critical exponents of the three- 
dimensional random-field Ising model are in agreement with the exponents 
of the pure Ising system in dimension 3 — 9 where 6 is the exponent that 
governs the hyperscaling violation in the random case. 
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The phase transition of Ising systems in a random quenched field is the subject of ongoing 
research Q. It is known that in the three-dimensional case, there is a low-temperature 
ordered phase. This was originally suggested by a heuristic argument due to Imry and Ma 
and has been shown rigorously @-[5|]. In the presence of random fields, these arguments 
show in fact that the lower critical dimension is d = 2 in the Ising case. The Hamiltonian 
of such systems is given by: 

H = -Y^SiS j - y £h i S i (1) 

where the Ising spins Si = ±1 are on the sites of a <i-dimensional cubic lattice and in- 
teract only between nearest-neighbors. The random fields hi have a quenched probability 
distribution hi = and hihj = h5ij where the overbar stands for average over the disorder. 
There is accumulating evidence for a second-order phase transition in three dimensions from 
Monte-Carlo studies real-space renormalization group calculations PHlTl as well as 

series expansions [O. 

The standard perturbative approach to critical phenomena leads to an upper critical 



dimension d = 6 for the random field problem and, below this dimension, the critical 
exponents are those of the pure Ising system in dimension d—2 within the e = d— 6 expansion. 
In fact this dimensional reduction holds to all orders in perturbation theory [0-[TH] in e. In 
view of the lower critical dimension d = 2, this cannot be the whole story. More general 



renormalization group studies |L7| have emphasized the role of a zero-temperature fixed- 
point ruling the physics at the transition point. This leads naturally to the presence of 
a new exponent 9 that governs hyperscaling violation: 2 — a = v(d — 9) where v is the 
correlation length exponent and a the specific heat exponent. This general scheme is in 
agreement with the droplet picture of the transition . In these theories, there is no 



obvious relationship between the random-field exponents as a function of the dimension and 
the exponents of the pure Ising system in dimension d — 9. Stated otherwise, the zero- 
temperature fixed point is a priori unrelated to the ordinary Ising fixed point. This implies 
also that there are a priori three independent exponents: 9 is a new exponent unrelated to 
v and 77, the exponent of the decay of spin correlations at the critical point. However, it 



has been suggested that dimensional reduction is still valid p2]j20|j2T|] in a generalized form, 
i.e. 9 = 2 — i] so that a two-exponent scaling is again valid. This means that the exponents 
of the random-field Ising model in dimension d would be those of the pure Ising model in 
dimension d — 9. 

In this Letter, we show that the best estimates for the critical exponents of the random- 
field Ising model in dimension d — 3 are in fact in agreement with the exponents of the 
pure Ising model analytically continued in dimension d' — 3 — 9. This analytic continuation 
is obtained through effective summation methods of the perturbative e = 4 — d! expansion 
series for the critical exponents of the pure Ising system. Such a procedure has been explored 
in the past and leads to very stable estimates fl22"| , j23| . 



The singular part of the free energy in a correlation volume F s i ng is governed by the 
correlation length: 

^sing CC e, (2) 
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where the spin correlation length diverges at the critical point £ oc \T — T c \ u . This definition 
leads to a modified hyperscaling relation: 

2-a = v{d-9). (3) 

Right at the critical temperature, the decay of the connected spin correlation function can 
be written: 

(S S X )-(S )(S X ) oc x-V-**} . (4) 
For the random field problem, the disconnected susceptibility has a different kind of scaling: 

(SoX^oc^. (5) 



This defines the two exponents 77 and 77. The general scaling theory of Bray and Moore f!7 
gives the following relation: 

9 = 2-r] + ri. (6) 



There is an exact inequality due to Schwartz and Soffer |24 



V < 2 V , (7) 

which in fact appears to be fulfilled as an equality 20|, 12| ,[7| . The scaling relations between 
exponents for the pure Ising model in dimension d! are obtained from those for the random- 
field case in dimension d through the relation: 

d' = d-e. (8) 

However, there is no direct evidence that the values of the exponents for the pure Ising model 
in dimension d' have any relationship with those for the random-field case in dimension d. 

To explore the relationship with the critical properties of the pure Ising model, we have 
used the estimates of the exponents as a function of the dimension d! coming from the 
e = 4 — d' expansion, up to order e 5 , of the pure system |23fl . To extrapolate the series, 
one needs to use powerful summation techniques that have been widely studied in the past 
P3| , |22|PB| . Let us explain here the principle of the procedure. If an exponent E is known 
from the e-expansion: 

E(e)=J2E h e k , (9) 

k 

one has to construct the so-called Borel transform B(et): 

POO 

E(e) = / dt exp{-t)t p B{et), (10) 



where p is an adjustable parameter. The series expansion for the Borel transform is then 
given by: 
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B{et) = Y,E k {tt) k /T{k + p + l) 1 (11) 

k 

where T is Euler's function. The coefficients are known to behave for large k as Ek ~ 
k\a k k h with some constants a, b. The advantage of the Borel transform is that this translates 
into a large k behaviour ~ a k k b ~ p for its coefficients in the series expansion. Thus the 
function B(et) is analytic at least in a circle: the closest singularity of B(eb) occurs at 
et = —1/a. Assuming that B(et) is analytic in the cut plane, one maps the cut plane onto 
a circle by use the mapping: 

^(T^-EW. (12) 

The Borel transform is now given through a convergent series in w, which leads to 

rod 

E(e) = J2B n ( dtexp(-t)t»[w(t)] n ). (13) 

Apparent convergence of such an expression is improved by varying p and by introducing 
other free parameters (see Ref. [p3~| ). This leads to stable estimates of the exponent E as 



a function of e. Errors in the extrapolation procedure can be estimated by varying the 
arbitrary parameters (like p) that enters the whole procedure. 

Extrapolating down from d! = 4 to d! = 2, one finds for example v = 0.99(4) for the 
correlation length exponent and 7 = 1.73(6) for the susceptibility exponent, in excellent 
agreement with the exact values from the Onsager solution v = 1 and 7 = 1.75. As 
a consequence, it is best to pin the exponents to the exact values in d! = 2 by writing 
E(e) = E(d' = 2) + (2 — e)E(e) and to perform then the above summation method on E. A 
global check of the method is provided by the comparison with other methods for d! = 3 : 
this scheme then leads for example to v — 0.6310(15) while the standard renormalization 
group result E5j from fixed dimension perturbation series is v = 0.6300(15), known to be in 



agreement with high temperature series and Monte-Carlo results. 

Reliable critical exponents 7, v, j3, 77 of the pure system have thus been obtained as a 
function of the dimension d! . In Figure 1, we plot the results for the susceptibility exponent 
7. The best value is the central solid line while the two extremal solid lines define the errors 
(there is thus a collapse at d! = 4 and d' = 2 of the three lines). Between d! — 4 and d' = 2 
the errors remain quite small and the value of the exponent does not vary very much. Below 
d! = 2, there is a rapid increase due to the proximity of the lower critical dimension of the 
pure problem. 

We then consider the random-field Ising model in dimension d = 3, for which reliable 
exponents are now available. The best Monte-Carlo results |7[ are obtained using a binary 
distribution of random-field. From the corresponding estimates of the optimal run with 
h/T = 0.35, one has t] = 0.56(3) and r\ = 1.00(6), and one then obtains from eq.(|8|) and 
eq.(P): d! = 1.44(10). The corresponding value for 7 is : 7 = 2.3(3). These Monte-Carlo 
results for the random-field Ising model in dimension d = 3 appear in Figure 1 as the solid 
box : its vertical size is fixed by the estimate of the exponent 7 with its error bar and its 
horizontal size is fixed by the range of d', which defines the solid box. 

In Figure 1, the large intersection between the solid box and the domain between the 
extremal solid lines, as well as the good intersection of the central solid line, clearly shows 
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FIG. 1. Critical exponent 7 of the magnetic susceptibility as a function of the dimension 
= 3 — 9. The meaning of the various symbols is given in the text. 



that 7 of the three dimensional random problem is in agreement with 7 of the pure system for 
d! = 3—8, and in particular close to 7 of the pure system for d! ~ 1.5 which corresponds to the 
simple guess |§ that random field fluctuations dominate the critical behaviour F sing ~ £ d / 2 : 
9 = 1.5 in three dimensions. We have also plotted in Figure 1 two less reliable other 
estimations: The dashed box is obtained in the same way as the solid one but with the results 
of the Monte-Carlo study of Ref. M using a Gaussian distribution of the random field, with 
a slower algorithm and a less extensive simulation than for the binary distribution. The two 
black dots are estimates from recent Migdal-Kadanoff renormalization-group studies |JT|D|,[TTJ] . 
Other results coming from earlier Monte-Carlo simulations |26| d' = 1.55(3) 7 = 1.7(2), from 
real-space renormalization-group calculations 7 = 1-9 — 2.2 or from high-temperature 
series expansion |12| 7 = 2.1(2) have not been drawn on Figure 1 for clarity, since they 
are well compatible with the quoted Monte-Carlo results. In Fig. 2, we plot the results 
for v with the same symbols. The dashed line is the estimate (up to d! « 1.6) coming 



from the near-planar interface model of Wallace and Zia |27| which has the same critical 
behaviour as the pure Ising model and can be expanded in powers of e' = d! — 1, giving 
v = e /_1 — 1/2 + e'/2 (We have plotted the only non-trivial Pade approximant). Note the 
agreement between this estimate and the e = 4 — d! expansion estimates. In Fig. 3, we 
plot the order parameter exponent f3. In this case, it is interesting to note that the Migdal- 
Kadanoff calculation ||To| , |i~T| shows a non-analytic vanishing of (3 near d = 2 and this also 
happens near d' — 1 in the droplet model of Bruce and Wallace p8| which is an extension of 
the near-planar interface model. Accordingly, our resummed results leads to an extremely 
small value (or even negative, the non-analytic behavior cannot be possibly reproduced by 
our approximants) of (3 for the corresponding dimension. 

Finally, in Fig. 4, we plot r/, again with the same symbols. In this CclS6 ; clS well as for /?, 
the Monte-Carlo results are in good agreement with the Migdal-Kadanoff values. 

The overall remarkable agreement clearly seen in our Figures between the results for the 
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FIG. 2. Critical exponent v of the correlation length as a function of the dimension d! . Same 
symbols as in Fig. 1 
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FIG. 4. Critical exponent rj of the spin correlations as a function of the dimension d! . Same 
symbols as in Fig. 1 

three dimensional random-field Ising model and those for the pure case for the corresponding 
dimension shows that the exponents of the random-field Ising model in three dimensions are 
in very good agreement with those of the pure Ising system in dimension d! = 3 — 9, and in 
particular very close to those of the pure system in dimension w 1,5 which corresponds to 
the a priori guess || 9 = 1.5 in three dimensions. 

We have thus obtained evidence for the generalized dimensional reduction for the random- 
field Ising model: the values of the exponents for the random model in dimension d are in 
fact those for the pure case in dimension d' = d — 9. 
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